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Abstract 


A propagator for the one dimensional time-dependent Schrodinger 
equation with an asymmetric rectangular potential is obtained using the 
multiple-scattering theory approach. It allows for the consideration of 
the reflection and transmission processes as the particle scattering at the 
potential jump (in contrast to the conventional wave-like picture) and 
for accounting for the nonclassical counterintuitive contribution of the 
backward-moving component of the wave packet attributed to the particle. 
This propagator completely resolves the corresponding time-dependent 
Schrodinger equation (defines the wave function ’tp{x,t)) and allows for 
considering the quantum mechanical effects of a particle reflection from 
the potential downward step/well and a particle tunneling through the 
potential barrier as a function of time. These results are related to fun¬ 
damental issues such as measuring time in quantum mechanics (tunneling 
time, time of arrival, dwell time). For imaginary time, which represents an 
inverse temperature {t —>■ —ih/ksT), the obtained propagator is equiva¬ 
lent to the density matrix for a particle that is in a heat bath and is subject 
to a rectangular potential. This density matrix provides information on 
the particles’ density in the different spatial areas relative to the potential 
location and on the quantum coherences of the different particle spatial 
states. If one shifts to imaginary time (t —>■ —it), the matrix element of 
the calculated propagator in the spatial basis provides a solution to the 
diffusion-like equation with a rectangular potential. The obtained exact 
results are presented as the integrals from elementary functions and thus 
allow for a numerical visualization of the probability density \tp{x,t)\'^, 
the density matrix and the solution of the diffusion-like equation. The 
results obtained may also be useful for spintronics applications due to 
the fact that the asymmetric (spin-dependent) rectangular potential can 
model the potential profile in layered magnetic nanostructures.. 

PACS numbers: 03.65.Nk, 03.65.Ta, 03.65.Xp 
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1 Introduction 

We start with the one-dimensional Schrodinger equation for a particle of mass 
m subject to potential V(x) 




where iT is a self-adjoint operator 


p2 


A solution to this equation can generally be presented as 


( 1 ) 

( 2 ) 

(3) 


where K(t) = exp{—iHt/h) is the propagator (Green’s function) for equation 
m in operator form and < x\K(t)\x' > is its matrix element in cc-representation. 
Thus, the knowledge of the propagator provides the complete solution to the 
equation dm at the given initial value 0). If the initial value is of the form 
'b(a;';0) = 6{x' — x"), the solution ([3]) reduces to the Green function matrix 
element 

'ip{x,x'';t) =< x\K{t)\x" > . (4) 
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Equation o with an imaginary time variable is also relevant to other phys¬ 
ical situations. If we make the substitutions t —>• —ihp (/3 = l/ksT) and 
'ip{x,x']—ihl3) —>■ p(a:, x';/3), Eq. (U) represents the matrix element p(a;, x';/3) =< 
x\ exp(—/3i7)|x' > of the density operator p{f5) = exp(—which satisfies the 
Bloch equation (in the x-representation) 


dp{x,x';P) 

dp 


Hp{x,x';P), 


(5) 


with initial condition p(x,x';/3 = 0) = S(x — x') and where the operator H ([2]) 
in (l5|) is applied only to the x variable of the density matrix. 

If we make the substitutions t —>■ —it, h —?> 2mD, V{x)/2mD —>■ V{x) 
and 'tp{x-,—it) —>• Q{x-,t), Eq. ([T]) represents the inhomogeneous diffusion-like 
equation (with the diffusion coefficient D) 


dQ{x;t) _ d'^Q{x;t) 

Wt “ (9x2 


V{x)Q{x-,t). 


( 6 ) 


The solution to Eq. ([6]) at the initial condition Q(x;0) = 6{x — xq) is given by 

O 

Q{x, xo; t) =< x\ exp(—i9t/2mi9)|xo >, (7) 

where H is defined by Q with h —>■ 2mD. 

We see that, in any case, the problem is to find a propagator of the type 
< x| exp(—ai7)|x' > with different a for the considered parabolic differential 
equations. 

A rectangular potential is the simplest one allowing for the study of some 
striking quantum mechanical effects, such as particle reflection from a potential 
step/well and transmission through a potential barrier. These phenomena are 
less surprising when we think of a wave being, e.g., reflected from a downward 
potential step, though they are more surprising from the particle point of view. 
They easily follow from the standard textbook stationary analysis, which re¬ 
duces to substituting a plane wave of energy E for the wave packet and solving 
the stationary Schrddinger equation. However, in this case, there are no real 
transport phenomena, i.e. in the absence of the energy dispersion [A.E = 0) 
the transmission time through or the time of arrival (TOA) to the potential 
jumps is indefinite (At ^ h/AE). It is interesting to verify the mentioned non- 
classical phenomena by considering the time-dependent picture of these pro¬ 
cesses in a realistic situation, when a particle, originally localized outside the 
potential well/barrier, moves towards the potential and experiences scattering 
at the potential jumps. In order to do this, the corresponding time-dependent 
Schrddinger equation needs to be solved. This problem is much more involved 
even in the one-dimensional case in comparison to the conventional stationary 
case. 

In particular, there is one striking and classically forbidden counterintu¬ 
itive (and often overlooked) effect even in the process of the simplest ID time- 
dependent scattering by the mentioned potentials. A wave packet representing 
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an ensemble of particles, confined initially (at t = to), say, somewhere to the 
region x < 0, consists of both positive and negative momentum components due 
to the fact that a particle cannot be completely localized at a; < 0 if the wave 
packet contains only p > 0 components. One would then expect that only par¬ 
ticles with positive momenta p may arrive at positive positions x > 0 at t > to- 
However, the wave packet’s negative momentum components (restricted to a 
half line in momentum space) are necessarily different from zero in the whole x 
space (— 00 -^ 00 ), representing the particles’ presence at x > 0 at initial moment 
of time to, and, therefore, may contribute, for example, to the distribution of 
the particles’ time of arrival (TOA) to x > 0 [TJ [5]. It is worth noting that 
the contribution of the backward-moving (negative momentum) components in 
the initial-value problem is in some sense equivalent to the contribution of the 
negative energy (evanescent) components in the source solution [1]. Thus, the 
correct treatment of some aspects of the kinetics of the wave packet (even in 
the ID case and even for a ’’free” motion) becomes a nontrivial problem and 
is closely related to the fundamental problem of measuring time in quantum 
mechanics, such as TOA, the dwell time, and tunneling time. 

In addition, the time-dependent aspects of reflection from and transmission 
through the potential step/barrier/well have recently acquired relevance not 
only in view of renewed interest in the fundamental problems of measuring time 
in quantum mechanics (see 0 ), but also due to important practical applica¬ 
tions in the newly emerged fields of nanoscience and nanotechnology. Rectan¬ 
gular (asymmetric spin-dependent) potential barriers/wells may often satisfac¬ 
torily approximate the one-dimensional potential profiles in layered magnetic 
nanostructures (with sharp interfaces). In such nanostructures, the giant mag¬ 
netoresistance (GMR) [4] and tunneling magnetoresistance (TMR) [5] effects 
occur. 

The calculation of the propagator < x\exp{—aH)\x' > is conveniently re¬ 
lated to the path-integral method (see, e.g. [6] and 0 )- The list of the exact 
solutions for this propagator is very short. For example, there is an exact so¬ 
lution for the spacetime propagator < x| exTp{—iHt/h)\x' > of the Schrodinger 
equation in the one-dimensional square barrier case obtained in [5] , but this so¬ 
lution is very complicated, implicit and not easy to analyze (see also [nmniin]). 

Recently, we have suggested a simple method for the calculation of the 
spacetime propagator [121 na [14], which exactly resolves the time-dependent 
Schrodinger equation with a rectangular potential in terms of integrals of ele¬ 
mentary functions. This method is an alternative to the commonly used path- 
integral approach to the mentioned problems and based on the energy integra¬ 
tion of the spectral density matrix (discontinuity of the energy-dependent Green 
function across the real energy axis). The energy-dependent Green function is 
then easily obtained for the step/barrier/well potentials with multiple-scattering 
theory (MST) using the effective energy-dependent potentials found in [T^ . 
which are responsible for reflection from and transmission through the potential 
step. These potentials, which are defined via the different particle velocities 
from both sides of the potential steps making up the step/barrier/well poten¬ 
tials, allow for the consideration of the reflection and transmission processes 
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as particle scattering at the potential jumps in contrast to the conventional 
wave-like picture. An important advantage of our approach is that the nega¬ 
tive energy (evanescent states) contribution to the propagator cancels out due 
to the natural decomposition of the propagator into forward- and backward- 
moving components. This is an essential result because accounting for both 
of these components (which should generally be done) often leads to a rather 
complicated consideration of the evanescent states with E < 0 (see [ID). 

In this paper, we provide an exact solution to Eq. (ED for real and imaginary 
times using our approach Eumn] to the calculation of the spacetime prop¬ 
agator for a general asymmetric rectangular potential. In Sec. 2, we outline 
our MST approach to the calculation of the propagator for the time-dependent 
Schrodinger equation and present its explicit form. In Sec. 3, we consider a 
system in a heat bath, as is the case, e.g., for electrons in nanostructures. The 
equilibrium system’s characteristics can then be calculated knowing its density 
matrix p{x,x';f3) =< x\exp{—f3H)\x' >. We present in this section an exact 
solution for the density matrix of a particle in an asymmetric (spin-dependent) 
one-dimensional rectangular potential and discuss its properties with the help of 
numerical evaluation of the corresponding integrals of elementary functions. In 
accordance with the above discussion of Eq. ED, the obtained solution for the 
spacetime propagator may be also used for finding the solution to the diffusion¬ 
like equation ED through the appropriate change of the equation parameters. 
This case is discussed in Sec. 4 and the summary of the results is given in Sec. 
5. 

2 Multiple scattering calculation of a spacetime 
propagator for the Schrodinger equation 

We start by considering a particle (electron) of mass m in the following gen¬ 
eral asymmetric one-dimensional rectangular potential of width d placed in the 
interval 0 < x < d 


V{x) = [e{x)-e{x-d)]U + 9{x-d)A, (8) 

where 0{x) is the Heaviside step function, and the potential parameters U and 
A can acquire positive as well as negative values (for A = U, V(x) reduces to 
the step potential). As an application, we can model by the potential ([HD a 
spin-dependent potential profile of a threelayer made of a nonmagnetic spacer 
(metallic or insulator) sandwiched between two magnetic (infinite) layers, and 
the asymmetricity (spin-dependence) of the potential ([HD is defined by the pa¬ 
rameter A. The particle wave vectors in different spatial areas (layers) are 
defined as 
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/ 2?77 

k%{E) = kiE),k{E) = < 0, 

kl{E) = k^^iE) = ku{E),ku{E) = ^‘^{E-U),0 <x<d, 

ki {E) = fcA (E),kA{E) = ^^{E-A),x>d. (9) 

In the case of three-dimensional sandwiches, and are the 

perpendicular-to-interface components of the wave vector k of a particle arriving 
at the interfaces (located at a: = 0 and x = d) from the right (>) or from the 
left (<). 

The wave function of a single particle moving in perturbing potential V(x) is 
given by Eq. 0 (see also i)- The propagator if (a;, a;'; t) =< x\exp{—iHt/H)\x' > 
is the probability amplitude for particle transition from the initial spacetime 
point (a:',0) to the hnal point {x,t) by means of all possible paths. It provides 
full information on the particle’s dynamics and resolves the corresponding time- 
dependent Schrodinger equation ([T]). According to [H], the time-dependent 
retarded propagator K{t) = 9{t — t') exp{—j-Ht) can be represented as 


K{t) = 9(t)— [ dEe-i^H --- 

27r 7 ^E-H + ie E - H - is 


),£—>■ - 1 - 0 , 


( 10 ) 


where H is the time-independent Hamiltonian of the system under considera¬ 
tion. Equation (1101) follows either from the contour integration in the complex 
plane or from the identity 


1 

E- H±ie 


= P T iT^5{E - H), 


( 11 ) 


where P is the symbol of the integral principal value. In the space representation 
(fTUl) reads 

OO 

K{x,x'-,t) = 9{t) j e~T^^*A{x,x';E)dE. (12) 

— OO 

Here A{x,x'-,E) is the spectral density matrix 


A{x, x';E) = ^ \G'^{x, x'\E) — G {x, x'; E)] , 

G^{x,x';E) =< x\——^ - —\x' >,G~{x,x';E) = [G~'~ {x', x; E)]* ,e —>■ -1-0, 

(13) 


determined by the matrix elements of the retarded G~^{E) and advanced G (E) 
energy-dependent operator Green functions = {E — H E ie)~^, which are 
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analytical in the upper and lower half-planes of the complex energy E, respec¬ 
tively. The propagator in the form of (IT^ is a useful tool for calculations within 
the multiple-scattering theory (MST) perturbation expansion if the Hamilto¬ 
nian can be split as H = Hq + Hi^ where Hq describes a free motion and Hi 
is the scattering potential. Note, that in this case one would not need rely 
on the standard (often cumbersome) matching procedure characteristic of the 
picture when a wave (representing a particle) is reflected from and transmitted 
through the potential ([5]). On the other hand, the introduction of the scattering 
potential Hi corresponds to the natural picture of the particle scattering at the 
potential jumps at a; = 0 and x = d. 

We showed in m that the Hamiltonian corresponding to the energy-conserving 
processes of scattering at potential steps can be presented as 


H = Ho + H,ix;E), 

H,{x-E)=Y,Ht{E)5{x-Xs). (14) 

S 


Here, Hi{x;E) describes the perturbation of the ’’free” particle motion (defined 
by Hq = — ^^t) localized at the potential steps with coordinates Xs (in the 
case of the potential ([5]), there are two potential steps at Xg = 0 and Xg = d) 


oh 

H!y{E) = -[vUE)-v^AE)]. 

H,<{E) = ^-^[v^^{E)-vUE)], 


H,><{E) 


2ihv?^{E)v^^{E) 


(15) 


where H^^^^-^{E) is the reflection (from the potential step at a: = a:s) potential 
amplitude, the index > (<) indicates the side on which the particle approaches 
the interface at x = Xg-. right (>) or left (<); Hf^^{E) is the transmission 
potential amplitude, and the velocities = hk^^^s^{E)/m, s G {0, d} 

(fc^(^)(i?) are given by (O). 

The perturbation expansion for the retarded Green function G+(a:, x'] E) in 
the case of the rectangular potential ([U, which can be effectively represented 
by the two-step effective scattering Hamiltonian (jTT)) . reads for different source 
(given by x') and destination (determined by x) areas of interest as 


G^{x, x']E) = G)}" {x, d; E)T^{E)Gq (0, x'; E),x' < 0,x > d, 

G+ix,x']E) = G^ix,0]E)T+{E)G^{d,x';E),x' >d,x<0, 

G+{x,x'-E) = G+{x,0;E)T'+{E)G^{0,x';E)+ G^ix,d;E)R'+{E)G^{0,x';E),x/ <0,0<x<d, 
G+(x,x';E) = G^(x,0;E)T'+(E)G^j-(0,x';E) + G;j-(x,0;E)H'+(E)G^j-(d,x';E),0< x^ <d,x<0, 
G+(x,x';E) = G^(x,x';E) +G;^(x,0;E)R+(E)G^(0,x';E),x' < 0,a: < 0, 

(16) 


7 





where the transmission and reflection matrices are 


_ Ti+^{E)G+{d,Q;E)Tl\{E) 

^ ’ D+{E) 

=T^+(i?)G+(d,0;£;)T'+(i?), 

^o+.z.^ , n+^{E)Gt(Q.d-E)Ti+{E)Gt{d,Q-E)Tl\{E) 

R {E)=T^ {E) + --, 

D+{E) = l-T^+{E)G+{d,0;E)T°+{E)G+{0,d-E). (17) 

The one-dimensional retarded Green function Gq{x,x'',E) corresponding to a 
free particle moving in constant potential V{x) =0 or V(x) = U{oi A) is (see, 

e.g. [12]) 

G^{x,x';E) = exp[ik{E)\x - x'\],V{x) =0, 

Goix,x;E) = ’^^(E’) exp[zfc^(A)(-G)k - x'\],V{x) = C/(or A), (18) 

where the wave numbers are determined by ([9]). The scattering (at the step lo¬ 
cated at a; = Xs) t-matrices are defined by the following perturbation expansion 


T^{E) = Ht{E) + Ht{E)Go{xs.x,;E)Ht{E) + ... 
HUE) 

l-Go{xs,Xs-,E)H^{E)^ 


(19) 


where EE’{E) and the interface Green function Go{xs,Xs',E) are defined differ¬ 
ently for reflection and transmission processes [12]: the step-localized effective 
potential is given by Eq. (ITSl) and the retarded Green functions at the interface 
for the considered reflection and transmission processes are, correspondingly, 


Xs', E) l/ihVy^^^{E) 

Goy<{xs,Xs-,E) = l/ih^v%{E)v^^{E) ( 20 ) 

in accordance with (fTSl) . 

From (ITSl) . (IT^ and (1201) . we have for the reflection T^'^^^(E) and transmis¬ 
sion T^~^(E) t-matrices, used in (fTTl) (s G {0,(i}) and corresponding to the re¬ 
tarded Green function and scattering at the interface located at x = Xg S {0, d}-. 

T;+iE) = ih./^t^, (21) 


where r^^^^(E) and t^(E) are the standard amplitudes for reflection to the right 










(left) of the potential step at x = Xg and transmission through this step 


XE) = 




-,r: 


US _ US 


t^{E) = 


2v^fcp^ 


(22) 


and the argument E is omitted for brevity. 

Now, using dsi), dni), (nzi), (HHi), dm) and (|2^ . we can obtain the Green 
function G'^{x,x'-,E) for the spatial domains considered in (IT6)) (see [16]) 


G+{x,x'-,E) = 




ih?\/kkA 
m 


x' < 0,x > d, 


G+{x,x'-,E) = , _ e-^^^t{E)E^^^^'-'^\x' > d,x < 0, 


ih?y/kkA 


G+{x,x']E) = 


m 


ih?\/kku - 


ikuX±f ( T^\ „ — ikx' I ^—ikuX„/{T^\„ — ikx' 


H'{E)e 


r'{E)e 


G+{x,x']E) = - 


ih?\/kku - 




G+{x,x'-,E) = - 


+r{E)e 


—ik{x-\-x') 


, x < 0, a;^ < 0, 


itflk . 

where the transmission and reflection amplitudes are defined as 

‘iy/kkXkue’-^'^’^ , 2y/M^{kA + ku) 

i(E) = -——- (E) = 


,x' < 0,0 < X < d, 
,x<0,0<x'<d, 
(23) 


d{E) 


d{E) 


, 2y/W^{ku - . {k-ku){kA + ku)-{k + ku){kA-ku)e^'^’^'^'^ 

^ yE) = - -77T^^ - ^tkE) = - 


d{E) 

d{E) = {k + fc„)(fcA + ku) - (k- fc„)(/cA - ku)e 


d{E) 


2ikud 


(24) 


Using the same approach, it is not difficult to obtain the Green function G'^ {x, x'] E) 
for other areas of arguments x and x'. 

In accordance with the obtained results for Green’s functions, we consider 
the situation when a particle, given originally by a wave packet localized to the 
left of the potential area, i.e. at x' < 0, moves towards the potential ([8]). We also 
choose A ^ 0, which corresponds to the case when, e.g., the spin-up electrons 
of the left magnetic layer [x' < 0) move through the nonmagnetic spacer to the 
right magnetic layer (x > d) aligned either in parallel (A = 0) or antiparallel 
(A > 0) to the left magnetic layer. At the same time, the amplitude U in the 
potential ([5]) may acquire both positive (barrier) and negative (well) values. 

From Eqs. (1231) we see that G+(x,x';A) = G+(x',x;A), and, therefore, 
the advanced Green function G~{x,x'-,E) = [G^ {x', x] E)]* = [G^(x, x'] E)]* 
(see, e.g. [E]). Thus, the transmission amplitude (fT^ is determined by the 
imaginary part of the Green function and can be written as 


E(x, x'-, t) 


TT 


J dEe-i^* linG+{x,x';E). 


— OO 


(25) 
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Formulas (1^ - (ESI) present the exact solution for the particle propagator in 
the presence of the potential ([8]). It should be kept in mind that the wave num¬ 
bers (HI) and, therefore, the quantities t{E), t'{E), r'{E), and r{E) in are 

0 oo 

different in the J dE and J dE energy integration areas: in the former case, 

— OO 0 

k{E) and k^iE) (A ^ 0) should be replaced with ik{E) and ik/\{E), where 
k{E) = \J—iLmEjk? (E < 0) and A:A(F';k||) = ■^2m{A. — E)/h?. At the same 
time, for energies E < 0, the wave number ku = iku, k^ = ^j2m{U — E^jh?^ 
for t/ > 0 (barrier), but for ?7 < 0 it is real, i.e. k^ = ^2m{E -(- \U\)/h^, if 
E > — \U\ and ku = iku, ku = ^J—2m{E + \U\)lh'^ li E < — \U\. It follows 
that the ’’free” Green function Gq {x, x'; E) = I is real in the energy 

interval (—oo-^O) and, therefore, does not contribute in this interval to the cor¬ 
responding ’’free” propagator Ko{x,x';t) defined by (1251) . It is also remarkable 
that for energies E < 0 the imaginary parts of the Green functions vanish in 
all spatial regions, as is seen from definitions dMl) and (IMl) (e.g., lmt{E) = 0 
and lmr{E) = 0 for E < 0). Therefore, the energy interval (—oo 0) does not 
contribute to the propagation of the particles through the potential well/barrier 
region. Thus, we have for t > 0 


K{x, x'; t) = 


I f dEe Ti 


— 4 Et 


■Xhj y/v{E) 


■ Re 


t{E)e 


(E)(x—d)^—ik{E)x' 




,x'<0,x>d, 


K{x,x'-,t) = — 
irh 


I 7dEe-i^* ( e-ik{E)x' U//^\^iku{E)x 

' ■ Re 






, x' < 0,0 < X < d. 


K{x,x'-t) = \ +r(A)e-*'=(^)("+"')],x' < 0,x < 0, 

Tin J v\E) 


irh J v{E) 

0 


(26) 


where the velocities v{E), Vu{E) and xa(A) are defined by ([9|) with the multi¬ 
plier h/m. 

It is easy to verify that the integration over E of the first term in the last 
line of (E51) results in the known formula for the space-time propagator for a 
freely moving particle 


Ko{x,x'-,t) 


/ \ll2 

<2rfS> 


l{x — x'Y 


2ht 


,x<0,x' < 0. 


(27) 


The obtained results (051) for the particle propagator completely resolve (by 
means of Eq. @) the time-dependent Schrodinger equation for a particle mov¬ 
ing under the influence of the rectangular potential ([8]). The form of this solution 
(integrals from the elementary functions) is convenient for numerical visualiza¬ 
tion. Further application of these results to the calculation of the TOA and 
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dwell time as well as of the probability density of finding a particle in different 
spatial areas as a function of time with account of the forward- and backward- 
moving components of the wave function and their interference can be found in 
our earlier papers [nnnmns]. 


3 Application to the density matrix 

The equilibrium non-normalized density operator (propagator in the temper¬ 
ature domain) p(/3) = exp(—/3iJ) can likewise be expressed in terms of the 
resolvent operator {E — H)~^ (see (ITU1) 1 

OO 

p(/3) = eM-m = f / dEe-P^[ ^ ^ 1 . )■ (28) 

zttJ E — H + le E — H — le 

— OO 

Particularly, in the coordinate representation we have for the density matrix 
(see (HI) 

OO 

p{x,x'-,f3)= j A{x,x'-, E)dE, (29) 

— OO 

where A{x,x'-,E) is given by ([13]). Thus, the density matrix p(x,x';/3) follows 
from the propagator m by the substitution t —>■ —ih(3 (/3 = l/kgT). From the 
properties (USD we see that the density matrix (1^^ is self-adjoint. The density 
operator (1^ satisfies the Bloch equation ([SJ. 

Thus, shifting to the imaginary ’’time” (t —>■ we obtain the exact 

density matrix p{x, x'; /3) in the various considered (relative to the potential ([8]) 
area) spatial regions, i.e. 


p{x, x'; /3) = K{x, x'; —ih/S), 


(30) 


where K{x,x'\t = —ihp) is given by (1261) . In particular, we obtain from (l27l) 
the known result for the ’’free” density matrix 


po(x, x\ /3) = Ko(x, x'; -ihjS) 


m 


'2T:h?-f5 


)V2exp[- 


m(x — x'Y 
2h'^l3 


(31) 


Using the same approach, it is not difficult to obtain the propagator p{x, x'; 13) 
for other (than in (1261) 1 areas of the arguments x and x'. Again, it is important 
to note that the negative energy half line (—oo^O), corresponding to the evanes¬ 
cent states does not contribute to the propagator ([29]). The diagonal element 
K{x, x] (3) {x = x' can be put only in the last line of (|26l) l defines the density of 
particles per unit length at the point a; < 0 to the left of the potential ([8|). The 
nondiagonal elements K{x,x'\ (3) of (1261) are related to the quantum mechani¬ 
cal interference effects. Particularly, they are responsible for particle tunneling 
through the barrier and also can be attributed to the phase correlation of the 
states I a: > and \x' >. 
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Equations (USD and (1501) provide an exact solution for the particle density 
matrix in the presence of the rectangular potential ([8]) in terms of integrals of 
elementary functions. It is convenient (e.g., for numerical visualization of the 
obtained results) to shift to dimensionless variables. As seen from ([5]), (jO]) and 
(jOOl) . there are the natural spatial scale d and the energy scale Ed = /2md? 

(the energy uncertainty due to particle localization within a potential range of 
width d). Then, the density matrix (1301) in the different spatial regions can be 
presented in the dimensionless variables as 


p{x,x'\P) = 


p{x,x'-,l3) = 


1 r dEe-P^ 

2TYd J E^/4: 

0 

oo ^ ' 

1 r dEe-P^ 

2 TYd J Ei/4: 


Re 


Re 


t{E)e 


\\/ E—A{x—l)^—i\/~Ex' 


,x'<0,x>l, 


{E - A)i/4 

t'(E)e^'^ _|_ 'E{E)e~^'^ e-Ux ^-iy/Ex 


(A - [/)i/4 


> , cc' < 0,0 < a: < 1, 


p(x, x'] P) = —exp[-(a; - x'f/4,p] + ^ 

2^d 


dEe-f^^ 

Ve 


Re[f(£;)e='^("+"')],J' < 0,2 < 0, 


(32) 


where 


t{E) 

t\E) 

r'{E) 


AEE^{E - Ap/^VE 

W) ’ 

2E^I^{E - Uy/^{VE-A + VE-U) 

W) ’ _ 

2E^I^{E - U)E^{^/e-U - Ve - 

W) 


r{E) = 


d{E) = 


{'/e- y/E - U){\/e -A + VE-U) - {y/E + \/E - U){v'e -A-^/E- 

^ ^ ^ ^ ^ ^ 

['/EE\lE-U){\jE-AE\lE-U) - {'/e - \J E - U){\/E - A - \J E - , 

(33) 


and E = E/Ed, U = U/Ed, A = A/Ed, P = Ed/keT = T/T, T = Ed/ke, 
X = x/d, x' = x'/d. 

We will visualize the results given by Eqs. (1551) and (1551) for several specific 
values of the relevant parameters. For an electron and the potential width d = 
10“"^CTO (Inm), the characteristic energy Ed ^ i- lQ~^ev and the characteristic 
temperature T = Ed/ks 3 • 10^AT. 

We will perform the numerical modeling of the density matrix (I32|) with the 
symmetric rectangular potential ([8]) when A = 0 (in this case the transition and 
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p(-2-2;10) 



Figure 1: The diagonal element p(—2,—2;10) as a function of the depth of the 


potential well 


U 


reflection amplitudes (l33l) simplify essentially). To secure a rapid convergence 
of the integrals in (15^ . we consider low enough temperatures, i.e. put /3 = 
10 (fcsT <C Ed). Figure 1 shows the diagonal element of the density matrix 
p{x,x;l3) (the last line in (l32ll l at a: = —2, i.e. the probability density to find 
a particle at this spatial point to the left of the barrier as a function of the 


potential well modulus 


U 


([/ = 0 


—300). We see that in this case the density 


matrix p{x,x',j5) exhibits a series of maximums and minimums. This can be 
explained by the formation of the resonance levels above the well if the condition 




U 


{n is integer, n = 1, 2,...) holds. With such a condition we have 


the reflection amplitude r(E) = 0 and the transmission amplitude t(E) = ±1. 
As at low temperatures (/3 = 10) the main contribution to the integral over E 
comes from the small (close to zero) energies, the positions of jumps at Fig. 1 


approximately follow the relation 


U 


= (n = 1, 2,...). 


The same diagonal element p(—2, —2; 10) as a function of the height of the 
potential barrier U = 0 ^ 300 behaves quite different from the case of the 
potential well and is shown in Fig. 2. One can see that the particle probability 
density at the given point to the left of the barrier x = —2 exhibits at firs^ a 
steep fall with the potential barrier growth and then it changes slowly with Lb 
We will evaluate the nondiagonal elements of the density matrix p(x,x';/3) 
for points at different sides of the well/barrier (the first line in (l32ll l as a function 
of the potential parameter U. Thus, we put x' = —10 (before the potential), 
X = 2 (beyond the potential) and /3 = 10 (as for Figs. (|S]) and ([5])). Figure 3 
exhibits the picks from the density matrix for the case of the potential well with 
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p(-2,-2;10) 



Figure 2: The same (as in Fig. 1) diagonal element of the density matrix as a 
function of the potential barrier height U. 


U = 0 —300 at the resonance values of 


U 


which correspond to the minima 


in Fig. 1. 

We see that in this case {U < 0) the density matrix nondiagonal elements 
can acquire both positive and negative values. Note that at [/ = 0 the density 
matrix reduces to the free density matrix dSU and therefore is positive (the 


seemingly negative value of the density matrix close to 17 = 0 in Fig. 3 is due 


to small resolution on the 


U 


-axis; calculation on the smaller scale near the 


point U = 0 shows that at 17 = 0 the density matrix is positive). Thus, Fig. 3 
demonstrates the jumps of the quantum coherence between the particle states 
before {x' < 0) and beyond {x > 0) the potential well at the resonance particle 


transmission through the potential well. At the values of 


U 


that do not satisfy 


the resonance condition this quantum coherence is small. 

The nondiagonal matrix element p(2, —10; 10) as a function of the potential 
barrier height (17 = 0 ^ 100) is shown in Fig. 4. We see that the quantum 
coherence between the states on the different sides of the barrier goes quickly 
enough to zero along with the barrier height. 


4 Diffusion-like equation 

As mentioned in the Introduction, the time-dependent Schrodinger equation o 
becomes equivalent to the parabolic diffusion-like equation ([6|) if one makes the 
substitutions t —>■ —it and h —^ 2mD, where H is a (diffusion) constant. Thus, 
we can immediately obtain from Eqs. (j26ll a solution to the diffusion equation 
([n|) for the initial condition Q{x-, 0) = 6{x — x'). In the dimensionless variables, 
this solution (a propagator) is given by Eqs. (l32l) . (p3l) with the following 


14 













p(2,-10;10) 



Figure 3: The nondiagonal element of the density 


function of the potential well depth 


U 


matrix p(2,—10;10) as a 


p(2,-10;10) 



Figure 4: The dependence of p(2, —10; 10) on the potential barrier height U. 
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Figure 5: The three-dimensional profile of the particles density Q{x, —3; t) to 
the right of the symmetric barrier. 


substitutions 


p{x, x'; (3) —>■ Q{x, X \ t), /3 —>■ t, i? —>■ i?, f/ —>■ 17, A —>■ A, 


t = t/tu, E = EjE d, U = U/Ed, A = A/Ed, 
tD = (f/D, Ed = 2mD^I(f, 


(34) 


where tD and Ed are obtained from td and Ed of the previous section by the 

substitution h —>■ 2mD. The introduced characteristic time Id, as it follows from 

the dehnition (IMl) , can be interpreted as the time needed for a particle to diffuse 

over the distance d (a potential (I5|) width) with the diffusion coefficient D. The 

2 

characteristic energy Ed = 2mD'^/d? = 2md^/t^ = 4^31^ jg proportional to 
the kinetic energy of a particle moving with the average velocity vd = d/tD- 
Therefore, as in the previous section, we can numerically model the solution 
Q(x, x'-, t) defined by Eqs. (15^ . (1551) (with the substitutions (1551) 1 in the different 
spatial points x = x/d and x' = x/d. 

Note that at V{x) >0 the solution to the diffusion equation ([6]) is positive, 
Q{x,x';t) > 0 (see [7]) and can be viewed as the density of particles in the 
point X at the moment of time t when the ’’diffusion with the holes” starts at 
the point x'. The latter term was introduced by Kac because in the points, 
where the potential ([5]) V{x) ^ 0 > 0); the particle can disappear. 

As an example, we have numerically modeled the density of particles Q{x, x';t) 
to the right of the symmetric barrier (A = 0) with [/ = 10 at different x, when 
the diffusion starts to the left of the barrier at x' < 0. The scaled time we chose 
t = 1 ^ 10 is sufficient to reach the spatial domain a; = 1 3 starting at x! = —3. 

The calculated three-dimensional profile of Q{x,x';t) is presented in Fig. 5 for 
the same (as earlier) width of the potential barrier d = 10~’^cm. 

One can see the nonmonotonic behavior of the density profile with time t 
for every fixed x, especially pronounced near the right barrier boundary (near 
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Figure 6: The spacetime density profile Q{x,—3-,t) for ’’free” diffusion of the 
particles. 


X = 1). This behavior, caused by the negative sources absorbing the particles 
(see Eq. ([6])) and distributed according to the function ([8|) with 17 > 0, A = 0, is 
quite different from the familiar ” free” diffusion in the absence of the potential 
{U = 0, A = 0) which is shown in Fig. 6 for the same parameters as in Fig. 5. 


5 Summary 

We have obtained the exact propagator < x\ exp(— ai7)|a;' > {H is the Hamilto¬ 
nian for a particle moving in the presence of the asymmetric rectangular poten¬ 
tial) resolving the parabolic-type partial differential equation. Having obtained 
the spacetime propagator for the one-dimensional time-dependent Schrodinger 
equation (a = it/h) with a rectangular well/barrier potential, we at the same 
time succeeded in finding a propagator for the Bloch equation [a = (3, j3 = 
l/kst) for the particle density matrix and for the diffusion-like equation {a = t) 
by shifting from real to imaginary time {t —>■ —ihfi and t —>■ —it, correspond¬ 
ingly). As an alternative to the conventional path integral approach to calculat¬ 
ing the propagators, we use the multiple-scattering theory for the calculation of 
the energy-dependent Green function (a resolvent operator in (llOp b The sug¬ 
gested approach is based on the possibility of introducing the effective potentials 
(see (HI and (I15p l which are responsible for reflection from and transmission 
through the potential jumps making up the rectangular potential ([5]). It pro¬ 
vides more of a non-classical picture of particle scattering at the considered 
potential as opposed to the conventional wave point of view. 

The solution for the time-dependent Schrodinger equation describes the re- 
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flection from and transmission through the asymmetric rectangular potential 
as a function of time and thus allows for considering the non-classical counter¬ 
intuitive effects of particle reflection from a potential well and transmission 
through a potential barrier in a real situation when a particle is moving towards 
the potential (|5|) and then experiencing a scattering at the potential. These re¬ 
sults are also relevant to the fundamental issues of measuring time in quantum 
mechanics such as the time-of-arrival (TOA), dwell time and tunneling time. 

The obtained density matrix p{x,x']j5) for a particle in a heat bath and 
under the influence of the potential (|S]) gives the probability density (diagonal 
matrix element) to find a particle in some spatial point and the quantum cor¬ 
relations (coherences) of different spatial states \x > and \x' > provided by the 
nondiagonal matrix elements. The results for the density matrix are numeri¬ 
cally visualized, which is enabled by the fact that they are expressed in terms 
of integrals of elementary functions. 

The results of the solution of the diffusion-like equation, which can be in¬ 
terpreted (for the case of a potential barrier, U > 0) as a diffusion with the 
negative sources distributed according to potential ([5]), also have been numeri¬ 
cally evaluated. The corresponding figures demonstrate the difference between 
this ” diffusion with the holes” and ” free” diffusion in the absence of the potential 
®. 

It is also worth mentioning that all obtained results are also relevant to 
the properties of electrons in nanostructures important for spintronics devices 
because the potential ® can be used for modeling potential profiles in such 
materials. 


References 

[1] A.D. Baute, I.L. Egusquiza, and J.G. Muga, J. Phys. A: Math. Theor. 34, 
42892001 (2001). 

[2] A.D. Baute, I.L. Egusquiza, and J.G. Muga, Int. J. Theor. Physics, 
Group Theory, and Nonlinear Optics 8, 1 (2002); (e-print arXiv.quant- 
ph/0007079). 

[3] J.G. Muga, R. Sala Mayato, and I.L. Egusquiza (ed). Time in Quantum 
Mechanics vol 1 (Lecture Notes in Physics vol 734) (Springer, Berlin, 2008). 

J.G. Muga, A. Ruschhaup, and A. del Gampo (ed). Time in Quantum 
Mechanics vol 2 (Lecture Notes in Physics vol 789) (Springer, Berlin, 2009). 

[4] M.N. Baibich, J.M. Broto, A. Fert, F. Nguyen Van Dau, F. Petroff, P. 
Etienne, G. Creuzet, A. Friederich, and J. Ghazelas, Phys. Rev. Lett. 61, 
2472 (1988). 

[5] R. Julliere, Phys. Lett. A 54, 225 (1975); 

P. LeGlair, J.S. Moodera, and R. Meservay, J. Appl. Phys. 76, 6546 (1994). 


18 


[6] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals 
(McGraw-Hill, New York, 1965). 

[7] M. Kac, Probability and Related Topics in Physical Sciencies (Lectures in 
Applied Mathematics vol 1) (Interscience Publishers, London, New York, 
1958). 

[8] A.O. Barut and I.H. Duru, Phys. Rev. A 38, 5906 (1988). 

[9] L.S. Schulman, Phys. Rev. Lett. 49, 599 (1982). 

[10] T.O. de Carvalho, Phys. Rev. A 47, 2562 (1993). 

[11] J.M. Yearsley, J. Phys. A: Math. Theor. 41, 285301 (2008). 

[12] V.F. Los and A.V. Los, J. Phys. A: Math. Theor. 43, 055304 (2010). 

[13] V.F. Los and A.V. Los, J. Phys. A: Math. Theor. 44, 215301 (2011). 

[14] V.F. Los and M.V. Los, J. Phys. A: Math. Theor. 45, 095302 (2012). 

[15] J.G. Muga and C.R. Leavens, Phys. Rep. 338, (2000). 

[16] V.F. Los and M.V. Los, Theor. Math. Phys. 177, 1704 (2013). 

[17] E.N. Economou, Green’s Eunctions in Quantum Physics (Springer, Berlin, 
1979). 


19 



